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Abstract 

A chain of one-dimensional Schrodinger operators connected by successive Dar- 
boux transformations is called the "Darboux chain" or "dressing chain". The pe- 
riodic dressing chain with period N has a control parameter a. If a / 0, the 
iV-periodic dressing chain may be thought of as a generalization of the fourth or 
fifth (depending on the parity of N) Painleve equations . The iV-periodic dressing 
chain has two different Lax representations due to Adler and to Noumi and Ya- 
mada. Adler's 2x2 Lax pair can be used to construct a transition matrix around 
the periodic lattice. One can thereby define an associated "spectral curve" and a 
set of Darboux coordinates called "spectral Darboux coordinates" . The equations 
of motion of the dressing chain can be converted to a Hamiltonian system in these 
Darboux coordinates. The symplectic structure of this Hamiltonian formalism turns 
out to be consistent with a Poisson structure previously studied by Veselov, Shabat, 
Noumi and Yamada. 
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1 Introduction 



A chain of one- dimensional Schrodinger operators L n , n G Z, connected by successive 
Darboux transformations is called a "Darboux chain" or "dressing chain" [SUES]- This 
is a kind of nonlinear lattice. A periodic dressing chain has a control parameter a. The 
nature of this system changes drastically whether a vanishes or not. If a = 0, the chain 
consists of finite-band operators EH Ej; the structure of those operators can be 
described by algebro-geometric methods 0112110]. If a ^ 0, the system is equivalent to 
the fourth and fifth Painleve equations and their higher order analogues [2H1I2]; the generic 
solutions are expected to be transcendental and beyond algebro-geometric methods. Thus 
the notion of periodic dressing chains invites us to an interesting world in which algebro- 
geometric functions and Painleve transcendents coexist and are connected by a continuous 
parameter. 

Periodic dressing chains have been studied from the point of view of Hamiltonian or 
Poisson structures as well. Veselov and Shabat [26J developed a bi-Hamiltonian formalism 
of periodic dressing chains. Noumi and Yamada jlH] took up substantially the same issue 
in the framework of their "higher order Painleve equations of type" (referred to as 
the Noumi- Yamada system in the following). Although discovered on a quite different 
ground, the Noumi- Yamada system of A^ type is actually equivalent to the periodic 
dressing chain with period N — i + 1. Noumi and Yamada define a Poisson structure on 
the phase space of this system, which coincides with one of the Poisson structures in the 
bi-Hamiltonian structure of Veselov and Shabat. 

On the other hand, a conceptually different Hamiltonian formalism has been devel- 
oped by Okamoto [22J for the Painleve equations and Garnier's generalizations [9.J. This 
Hamiltonian formalism is based on the notion of isomonodromic deformations of an or- 
dinary differential equation. Harnad and Wisse fOl E] pointed out that the notion of 
"spectral Darboux coordinates" pQ lies in the heart of this kind of Hamiltonian formal- 
ism of isomonodromic deformations. This reveals a close relationship to "separation of 
variables" of integrable systems pHl EH] . Since periodic dressing chains include the fourth 
and fifth Painleve equations, one will naturally ask whether periodic dressing chains have 
a similar Hamiltonian formalism. This is the main problem that we address in this paper. 

For this purpose, we need an isomonodromic Lax formalism of periodic dressing chains. 
Two options can be found in the literature. One is a Lax pair of 2 x 2 matrices used by 
Veselov and Shabat [2E] for the case of a = and modified by Adler [H] to fit into the case 
with an arbitrary value of a. Strictly speaking, Adler 's Lax pair is related to deformations 
of a difference (rather than differential) equation. Therefore this is not an isomonodromic 
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system in the usual sense. Another option is a Lax pair of N x N matrices presented by 
Noumi and Yamada 19j. This is an isomonodromic system in the usual sense. These two 
Lax pairs turn out to be equivalent (or dual) and connected with each other by a Mellin 
transformation. 

We use Adler's Lax pair and construct a 2 x 2 transition matrix T(A) around the 
periodic chain. This matrix satisfies a Lax equation, which inherits the unusual nature 
of Adler's Lax pair. Nevertheless, we can construct a hyperelliptic "spectral curve" T 
thereof, and convert the Lax equation to a dynamical system of a finite number of moving 
points (Xj,Zj) 6 r, j = 1, . . . , g. Another complicated aspect of this system is that the 
spectral curve itself is dynamical in the case where cr ^ 0, but this is rather a universal 
characteristic of isomonodromic systems. The coordinates Xj, Zj of the moving points are 
nothing but "spectral Darboux coordinates" mentioned above. The equations of motion 
of the moving points becomes a Hamiltonian system (which we call, for convenience, the 
spectral Hamiltonian system) in these new coordinates, and conversely, the Lax equation 
can be reconstructed from this Hamiltonian system. 

It is not obvious from the construction that the passage from the periodic dressing 
chain to the spectral Hamiltonian system is invertible. The problem is to reconstruct 
the phase space coordinates of the periodic dressing chain or, equivalently, of the Noumi- 
Yamada system from the spectral Darboux coordinates. We solve this inverse problem 
and find the existence of a locally invertible map that connects the spectral Hamiltonian 
system with the Noumi- Yamada system directly. As a byproduct, we can show that 
this map is a Poisson map connecting two apparently different odd-dimensional Poisson 
structures. 

This paper is organized as follows. Sections 2 and 3 review the equations of motion 
of dressing chains, the Noumi- Yamada systems, and Lax pairs of these systems. Section 
4 deals with the transition matrix and the associated spectral curve. Sections 5 and 6 
present the main results of this paper. Section 5 is concerned with the Hamiltonian system 
in spectral Darboux coordinates. Section 6 is devoted to the inverse problem and Poisson 
structures. Section 7 is for concluding remarks. Some technical details are collected in 
Appendices. 
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2 Dressing chain 
2.1 Definition 

A dressing chain consists of one-dimensional Schrodinger operators of the factorized form 

L„ = (d x - v n )(d x + v n ), d x = d/dx, (2.1) 

that are linked with the neighbors by the Darboux transformations 

(d x - v n+1 )(d x + v n+1 ) = (d x + v n )(d x - v n ) + a n (2.2) 

with parameters a n . This is a kind of nonlinear lattice governed by the differential equa- 
tions 

v n + v n+ i = v n+1 -v 2 n + a n , (2.3) 
where the dot stands for the x-derivative: 

Vn d x V n - 

We are mostly interested in periodic dressing chains satisfying the periodicity condition 

v n+N = v n , a n+N = a n . (2.4) 

Since the case of iV = 2 is no very interesting, the subsequent consideration is focussed 
on the case where N >3. The auxiliary variable 

N ^ N 

v = ^2v n = -^2(v n + V n+ l) (2.5) 

n=l A n=l 

satisfies the very simple equation of motion 



v = g , (2.6) 



where a is the constant defined by 
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a = J2 a n- (2-7) 



n=l 



Thus v is a linear function of the form (a/2)x + constant. This means that v is not a 
dynamical variable in the usual sense. 



The parameter a is a control parameter of the periodic dressing chain. To see how 
the situation depends on this parameter, we note that the iV-th order operator 

M n = (d x + Vn+N-x) ■ ■ ■ (d x + O (2.8) 

satisfies the equation 

[L nj M n ]=aM n . (2.9) 

If a = 0, L n and M n are commuting differential operators in the sense of Burchnall and 
Chaundy [4 . In particular, the case of odd iV's falls into the classical theory of finite-band 
operators. The case of even iV's is more involved, because the orders of the operators are 
not co-prime. If a ^ 0, the operator equation resembles the so called "string equations" 
of non-critical string theory and two-dimensional quantum gravity ^H] . The usual string 
equations take the form 

[L,M]=a (2.10) 

(scaling limit of Hermitian matrix models) or 

[L,M]=aL (2.11) 

(scaling limit of unitary matrix models). In particular, if L is a second order Schrodinger 
operator L = d 2 x — u, the string equations of the first type give the first Painleve equation 
and its higher order analogues. The string equation of the second type are related to the 
second Painleve equation. Similarly, the periodic dressing chains for N = 3 and N = 4 
are equivalent to the fourth and fifth Painleve equations, as pointed out by Veselov and 
Shabat and Adler 0. 

2.2 Relation to Noumi-Yamada systems 

The periodic chains have another expression in terms of the new dependent variables 

fn = V n+ i+V n . (2.12) 

As we shall see below, this gives a change of variables to the Noumi-Yamada system [18^ 
2.2.1 If N is odd 

If N = 2g + 1, ()2.12j) can be solved for v n as 

Vn = gC/tl - fn+1 H h fn+2g)- (2.13) 
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In this case, t> n 's and / n 's are thus connected by an invertible linear map. Rewriting ()2.3|) 

as 

fn = fn(v n +l - V n ) + a n 

and inserting this expression of v n therein, one obtains the closed differential equations 

/ 9 9 \ 



fn = fn[Yl fn+2k-l ~ f n + 2k ) + a " ( 2 - 14 ) 

Vfc=i k=i ) 

for / n 's. This is exactly the Noumi-Yamada system of A^j type. 



2.2.2 If N is even 

If N = 2g + 2, (|2.12j) cannot be solved for v n 's uniquely. This ambiguity, however, turns 
out to be removed by the existence of an extra constraint. This constraint is a consequence 
of the obvious identity 

29+2 

E (-!)"(«» + «n+i) = (2.15) 

n=l 

and the equations of motion of the dressing chain. Because of this extra constraint, v n 's 
and fnS are connected by a birational map; see Appendix B for details. Thus one can, in 
principle, rewrite the equations of motion of the dressing chain into a system of differential 
equations for / n 's. We can, however, circumvent messy calculations that will be inevitable 
if we do the change of variables naively, as follows. 

The clue is to calculate the derivative of vf n rather than /„ itself. Since 

9+1 

V = ^2 fn+2k-l, 
k=l 

the derivative of vf n can be written as 

9+1 



V 



fn) — ^,{fk+2k-lfn + fn+2k-lfn 



k=l 



We can now use ()2.3|) to rewrite each term of the sum on the right hand side as 

fn+2k-lfn + fn+2k-lfn 

— (fn+2k-l( v n+2k ~ V n+ 2k-l) + «n+2fc-l )fn fn+ 

= /n+2fc-l/n(/n+2fc-l — 2fn+2k-2 H ' + 2fn+l ~ fn) 

+ &n+2k-lfn + (*nfn+2k-l- 
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Suming over n and doing some algebra, we eventually obtain the differential equations 

( v fn) — fn I £ fn+2j-lfn+2k ~ £ fn+2j fn+2k+l 
\l<j<k<g l<j<j<9 

9+1 

+ £ a n+2 k-ifn + a n v (2. 16) 

fc=i 

or, equivalently, 

^/n = /n I X! fn+2j-lfn+2k ~ £ fn+2jfn+2l,+l 
\l<j<k<g l<j<fc<9 

+ I £ a n+2fc-l - - I /» + a n V. (2.17) 

If f is normalized as v = ax/2, this is nothing but the Noumi-Yamada system of type 
^20+1 ■ Note that the variables f n , being defined by ()2.12|) . also satisfy the constraint 

9+1 9+1 

£/tt-i = £/2*- (2-18) 

fc=i fe=i 

3 Lax pairs 

The dressing chains have two different Lax pairs, one due to Veselov and Shabat |2E] and 
Adler [3j, and the other due to Noumi and Yamada We point out that these two 
Lax pairs are connected by a Mellin transformation. 

3.1 2x2 Lax pair 

The Lax pair of Veselov and Shabat is valid only for the case of a = 0. Adler proposed a 
slightly modified Lax pair in which this assumption is removed. 
Let u n denote the function 

u n = v 2 n -v n . (3.1) 
This is the "potential" of the Schrodinger perator L n : 

L n = d 2 x - u n . (3.2) 
Adler's Lax pair consists of the 2x2 matrices 

V n {\)=\ *" 2 1 ) (3.3) 




and 

u n (\)= k ; ; i , (3.4) 

which satisfy the differential equation 

V n (X) = U n+1 (X + a n )V n (X) - V n (X)U n (X). (3.5) 

This "Lax" equation resembles the Lax equation of many integrable nonlinear lattices; an 
unusual feature is that it contains a shift of the spectral parameter. 
An associated auxiliary linear system takes the form 

$ n (A) = C/ n (A)$ n (A), (3.6) 
$ n+1 (A + « n ) + K(A)$„(A) = 0. (3.7) 

The Lax equation ()3.5|) can be reconstructed from these linear equations as the consistency 
(integrability) condition. One can rewrite this matrix linear system to a scalar form. 
Firstly, (|3.fij) implies that $ n (A) can be expressed as 



$„(A) 



4>n(X) 

MX) 



and that the scalar function <p n {X) satisfies the Schrodinger equation 

L n cf) n (X) = 0„(A) - u n (j) n (X) = A0 n (A). (3.8) 
Secondly, ()3.7|) yields the the scalar difference-differential equation 

n+ i(A + an) + MX) + v n M^) = 0- (3.9) 

3.2 Another Lax pair 

Another Lax pair emerges from the auxiliary linear system 

azd z ip n (z) + ip n+2 {z) + f n ip n+ i{z) + e n tp n (z) = 0, (3.10) 

1p n (z) + 4> n +l(z) + V n 1p n {z) = 0. (3.11) 

Since (|3.11j) implies that 

= (d x + v n+ i)(d x + v n )i) n (z) - (v n + v n+1 )(d x + v n )if) n + e n i) n (z) 
= (L n + e n )ip n (z), 



8 



one can rewrite (j3.10|) as 



azd z ip n (z) + (L n + e n )ip n (z) = 0. (3.12) 
The consistency condition of these linear equations yields the three equations 

fn+i ~ fn = v n+2 - v n , (3.13) 

fn = fn{Vn+l ~ V n ) + C n - e n+1 , (3-14) 

e n = 0. (3.15) 
These equations are retained under gauge transformations 

^n^^Pne 1 , V n ^V n -j, fn^fn, e„ -> C n . (3.16) 

The equations of motion ()2.3)1 of the infinite dressing chain can be recovered by a special 
"gauge fixing" in which ()2.12j) is satisfied. Let us rewrite the first equation ()3.13j) of the 
consistency condition as 

fn+l — V n +2 ~ v n+l — fn ~ ~ V n . 

Both hand sides of this equation is thus independent of n, thereby become a function h = 
h(x) of x only. This function transforms as h — > h + 27 under the gauge transformation, 
and can be eliminated (or "gauged away") by a suitable choice of the function 7. One 
can thus achieve a special gauge that satisfies ()2.12j) . In this gauge, the second equation 
()3.14|) of the consistency condition becomes the equation ()2.3|) of the dressing chain upon 
identifying the parameters a n as 

a n = e n -e n+1 . (3.17) 

Note that a is so far a free parameter; the relation to a n 's emerges after imposing the 
iV-periodicity condition (see below.) 

This auxiliary linear system reduces to that of Noumi and Yamada ^Hj for higher 
order Painleve equations by imposing the periodicity conditions 

Vn+N = Vn, fn+N = fn, &n+N = C n — Oi. (3.18) 

The last condition implies the relation 

N 

a = J2(e n 

n=l 

9 



N 

- e n+ i) = a n 

n=l 



as expected. We further assume the quasi-periodicity 

4>n+N(z) = Zlf) n (z) 



(3.19) 



of ip n (z), z being interpreted to be the Bloch-Floquet multiplier. The auxiliary linear 
system of the infinite chain thereby turns into the N x N matrix equations 



azd z V(z) + £(z)ty{z) = 0, + M(z)^(z) = 
for the column vector 

= \^ x {z) ••• M< 
where C(z) and M.(z) denotes the N x N matrices 

( ei fx 1 



(3.20) 



C(z) 



M(z) 



V fNZ Z 
I Vi 1 



V z 



(3.21) 



6jV— 1 /aT-1 

\ 



1 



(3.22) 



This Lax pair may be thought of as a kind of self-similarity reduction of the Drinfeld- 
Sokolov systems [Sj. 

3.3 Mellin transformation connecting two Lax pairs 

The two auxiliary linear systems turn out to be connected by the Mellin transformation 

Mz) = J ^" AA >«(A - e n )d\ (3.23) 
provided that (|3.11|) can be derived from ()3.9|) as 

ip n +i{z) = J z~ x/a (j) n+l (\ - e n + a n )d\ 

= - J z- x ' a {d x + v n )(t> n {\-e n )d\ 
= -(d x + v n )ip n (z). 
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Similarly, (j3.12j) is a consequence of the identity 



azd z ip n (z) 




and the equation 



A</> n (A - e n ) 



(A - e n )(p n (X - e n ) + e n <p n (X - e n ) 
(L n + e n )0 n (A - e n ) 



that can be derived from (J3.8)) . 

Note that the Mellin transformation is a heuristic relation rather than a rigorous 
one. For instance, we have not specified the meaning of the integral. Nevertheless, this 
correspondence is very useful for understanding the origin of the two different Lax pairs. 
Also note that the Mellin transformation can be converted to a Laplace transformation 
by changing the variable from z to ( = log z. 

4 Transition matrix and spectral curve 

We construct the transition matrix T(A) of periodic dressing chains. This matrix satisfies 
a Lax equation, which is not isospectral if a 7^ 0. The spectral curve V is accordingly 
"time-dependent", i.e., deforms as x varies. 

4.1 Transition matrix 

Following the usual prescription, we can consider an analogue of the transition matrices 
around the periodic lattice. The transition matrix T(A) that connecting $i with §n+i is 
defined by the linear relation 



Note that A is also shifted as A — > A + a. More explicitly, T(A) can be written as 



$ w+1 (A + «) = (-l) JV T(A)$ 1 (A). 



(4.1) 



T(A) = V N (X + ,%-!) ■ • • V 2 (X + A)K(A) 



(4.2) 



where 



Pn = «i + «2 H ha n , (3 = 



(4.3) 



and satisfies the Lax equation 



T(A) = U X {X + a)T(X) - T(X)U 1 (X). 



(4.4) 
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Note that this Lax equation is not isospectral if a ^ 0. In terms of the matrix elements, 

T(A) = 



A(X) B(X) 
C(X) D(X) 



the Lax equation of T(A) reads 

A(X) = C(X)-B(X)(X + Ul ), (4.5) 

B(X) = D(X)-A(X), (4.6) 

C(X) = (X + a + Ul )A(X)-D(X)(X + Ul ), (4.7) 

D(X) = (X + a + Ul )B(X)-C(X). (4.8) 

4.2 Equation of spectral curve 

The affine part of the spectral curve T is defined by the equation 

det(zl - T(A)) = z 2 - P{X)z + Q(X) = (4.9) 
on the (A, z) plane. -P(A) and Q(X) are the spectral invariants 

P(A) = TrT(A), Q(A) = detT(A) (4.10) 
of T(A). Since T(A) is a product of 2 x 2 matrices, Q(X) can be readily determined: 

N-l 

Q(X) = (-l) N XY[(X + p n ). (4.11) 

71=1 

To see what P(X) looks like, one has to specify the matrix elements of T(A) in more 
detail. The following can be confirmed rather easily by induction on iV: 

1. If N — 2g + 1, the matrix elements of T(A) are polynomials of the form 

A(X) = a X 9 + aiX 9 ' 1 H h a g , 

B(X) = b X° + b 1 X^ + --- + b g , 
C(A) = c X^ + Cl X 9 + ■■■ + %+!, 
D(X) = d X 9 + diX 3 - 1 + • • • + d g , 

where 

bo = Co = 1, a = do = v. (4-13) 
Consequently, -P(A) is a polynomial of the form 

P(X) = I X 9 + hX 9 ' 1 + ■ ■ ■ + I g . (4.14) 
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2. If N = 2g + 2, the matrix elements of T(A) are polynomials of the form 

A(X) = X° +1 + a \9 + ■ ■ ■ + a g , 
B(X) = b \v + b 1 \°~ 1 + --- + b g , 
C(X) = c \° +1 + Cl \ 9 + ■ ■ ■ + c g+1 , 
D(X) = X 9+1 + d X 9 + dxX 9 - 1 + • • • + d g , 

where 

b = c = v. (4.16) 
Consequently, P(X) is a polynomial of the form 

P(A) = 2X 9+1 + J A 9 + h\ 9 ~ l + ■ • • + J - (4.17) 

T is a hyperelliptic curve of genus g. The defining equation can be transformed to the 
normal form 

y 2 = P{Xf-AQ{X) (4.18) 

by changing the variables as z = (-P(A) + y)/2. If a = 0, the coefficients I n of -P(A) are 
constants of motion; if a ^ 0, they depend on x. This is obvious from the equation 

P(A) = aB(X), (4.19) 

which can be derived from ()4.5|) and ()4.8)1 . 
A consequence of (I4.19J1 is the equation 

J = ab . (4.20) 

Since b is a constant (for odd iV's) or a simple linear function of x (for even A^'s), this 
equation shows that J Q , too, is a simple (at most quadratic) function of x even if a ^ 0. 
We shall present a more explicit expression of Iq below. 

One can now find an explicit expression of the coefficients a® and do for iV = 2g + 2. 
This is a consequence of ()4.6| : this equation implies that 

b = d - a . (4.21) 

Since d + ao = lo and 6o = v = a/2, this immediately leads to the expression 

«o = ^(/o " 6d) - \{h ~ f ), do = i(/o + 6d) = ^(/o + f ) (4-22) 

of ao and do in terms of Jo- 

This expression implies that a and do, like bo, Co and Jo, are a simple (at most 
quadratic) function of x (in particular, a constant if a = 0). We shall interpret these 
quantities as Casimir functions of a Poisson structure. 
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4.3 Structure of P(A) in more detail 

We here review an explicit formula of -P(A) presented by Veselov and Shabat 
By definition, P(X) is given by the trace 

P(A) = Tr V N (X + (3 N -i) ■ ■ ■ V 2 (X + /W(A). (4.23) 

We now split each matrix V n (X + in the trace into two rank-one matrices: 

w+a-o= („!)("» o+^Jo °)- (4 - 24) 

For convenience, let us call the two factors on the right hand side as type I and type II 
factors, respectively. 

The trace is now expanded into a sum of 2 N terms, each term being the trace of a 
product of rank-one factors as above. The trace of a product of such matrices factorizes 
in itself. For instance, if all the rank-one factors in the trace is of type I, 



Tr ( \{v N l)---\ \[v 1 l) = (v N + vn-!) ■ ■ ■ (v 2 + fi)Oi + v N ) 

= /1/2 • - ■ /iV- 

The expansion of P(A) starts with this term, and proceeds to higher order terms including 
type II factors. General terms obey the following rules. 

1. A higher order term in this expansion vanishes if it includes a pair of neighboring 
factors of type II. This is due to the obvious identity 

a + p n -i J y A + l3 n -2 J 1 \ 

2. After factorization, non- vanishing traces take place in the following three forms: 



Tr ( v n+ i l)( 

V ' \ V n 



Tr( v n+1 1 ) ( A+ °^ ^ J = A + /3 n _ x , 



1. 



In particular, v n always appears along with v n ±! in the linear combinations v r , 
v n +i = fn and v n + u n _i = / n _i. 
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This expansion can be reorganized into the beautiful formula 

N-l / 02 \ 

PW= II 1 + (A + /3 TO ) )(fif2---f N ) (4.25) 

n=0 V OJnOJn+lJ 

due to Veselov and Shabat [2E]. Here it is understood that fo = and f3 = 0. One 
can read off some of the coefficients of -P(A) explicitly. In particular, J turns out to be a 
polynomial in v with constant coefficients: 

1. UN = 2g + l, 

lo = fi + h + ■ ■ ■ + h 9 +i = 2v. (4.26) 

2. UN = 2g + 2, 

g+l 2g+l 2g+l 

h = E hi-ihk + E Pn = v 2 + E (3 n . (4.27) 

j,k=l n=l n=l 

5 Spectral Darboux coordinates and Hamiltonian sys- 
tem 

The Lax equation of the transition matrix can be converted to a dynamical system of 
a finite number of points on the spectral curve. The coordinates of these points are 
"spectral Darboux coordinates" in the terminology of the Montreal group mElJE]. ^ ne 
dynamical system turns out to take a Hamiltonian form in these Darboux coordinates. 

5.1 Spectrl Darboux coordinates 

Let Xj, j = 1, ■ ■ • , g, denote the zeroes of -B(A): 

B(\) = b f[(X-X J ). (5.1) 

i=i 

We consider the generic case where the zeroes A-,- are mutually distinct, and derive a 
system of differential equations for these zeroes. 

Differentiating the identity B(Xj) = by x yields the identity 

B(Xj) + B'(Xj)Xj = 0, 

where B'(X) stands for the A-derivative: 

B'(X) = d x B(X). 
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Solving this identity for Xj and recalling ()4.6|) . we have 

In order to rewrite the numerator of this equation in a more familiar form, let us notice 
the identity 

(A(X) - D(X)) 2 = P(X) 2 - AA(X)D(X). 
Upon setting A = Xj, this identity becomes 

(A(Xj) - D(Xj)f = P(Xjf - 4Q(Xj) 

because 

A(Xj)D(Xj) = A(Xj)D(Xj) - B(Xj)C(Xj) = Q(Xj). (5.3) 
Thus A(Xj) — D(Xj) turns out to be expressed as 



A(Xj) - D(Xj) = y/P{\ s )* - 4Q(Xj). (5.4) 
This leads to the differential equation 



t 'W-4q(A,) .... 

- BW (5 - 5> 

for the zeroes of -B(A). 

The last equation resembles the so called "Dubrovin equation" in the theory of finite- 
band integration [2H EH E]- I n fact, if a = 0, this is indeed a variant of Dubrovin 
equation that can be solved by the usual algebro-geometric method. If a ^ 0, the situation 
drastically changes. Namely, the coefficients of -P(A), as well as bo, are no longer constant, 
r itself is accordingly a dynamical object, for which (|4.19J) plays the role of equations of 
motion. 

The roots Ai, . . . , X g of B(X) comprise a half of the spectral Darboux coordinates. The 
other half defined by 

z 3 = A(Xj). (5.6) 

The g-tuple of points (Xj, Zj) of the (A, z) plane sit on the spectral curve Y, namely, satisfy 
the algebraic relations 

z 2 -P(X J )z 1 +Q(X j ) = 0, (5.7) 
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as one can readily see from the triangularity 

T(A,) = 



A(\j) 



of T(Xj). We can prove that these moving points on T satisfy the following differential 
equations. 

Proposition 1 Xj and zj satisfy the differential equations 

• _ Zj - Q{X J )zj 1 



P'(X,)z, - Q\X 3 ) 
*i ~ s\X~, • (5.9) 



Proof. (15.8)) is nothing but (J5.2j) in disguise; the two terms in the denominator of (J5.2j) 
can be rewritten as 

A(X j )=z ji D{X 3 ) = ^ ) =Q{X j )zj\ 

To derive (|5.9|l . we differentiate (|5.7jl by x. This yields 

(2 Zj - P\X 3 )) Zj + (-P'(Xj)zj + Q'(Xj))Xj - P(Xj)zj = 0. 

Solving this equation for Zj and noting the algebraic relation P{Xj) = that holds by 
(|4.19|) . we find the equation 

. P'(Xj)zj - Q'{Xj) • 
2i " 2 Zj - PiXj) 

On the other hand, (|5.2jl can also be written as 

• 2A(Aj) - PjXj) _ 2z 3 - P(\j) 
B'iXj) B'(Xj) ■ 

Eliminating Xj from these two equations, we obtain (J5.9j) . □ 



5.2 Deriving Hamiltonian system 

We now seek to convert the differential equations of Xj and Zj to a Hamiltonian system 
(which we call the "spectral Hamiltonian system"). Since the situation is parallel to the 
case of the periodic Toda chains jH], we first review that case, then turn to the case of 
periodic dressing chains. 
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5.2.1 Periodic Toda chains 

The description of periodic Toda chains, too, starts from a 2 x 2 Lax pair and an associated 
transition matrix. The spectral curve of an iV-periodic chain is a hyperelliptic curve of 
genus g = N — 1 defined by the equation 

z 2 - P(X)z + Q = 0, (5.10) 

where P(X) and Q are the trace and the determinant of the transition matrix. Whereas 
P(X) is a polynomial of the form 

N 

P(\) = \ N + Y,Ii\ N - e , (5.11) 
e=i 

Q is now a constant. The next-to-leading term I x vanishes in the center-of-mass frame, 
and the other coefficients I2, • ■ ■ , In are constants of motion; I2 can be identified with 
the Toda Hamiltonian H. The first two matrix elements A(X) and B(X) of the transition 
matrix determine a set of moving points (Ai, Z\), ■ • ■ , (Ajv-i, zn-x) on the spectral curve as 
B>(Xj) = and zj = A(Xj). Consequently, the coordinates of these moving points satisfy 
the algebraic equations 

z) - P{Xj)zj + Q = 0. (5.12) 
Using the inversion formula (jA.5|) in Appendix A, one can solve these equations for In as 

h ~ ~ g *m ox/ (5 ' 13) 

where bn are the coefficients of B(X), 

B(X) = X N - 1 + b 1 X N - 2 + --- + b N _ 1 , 

the leading coefficient being equal to 1 in the usual setup. Since b\ = — Ai — ■ • • — Ajv-i, 
we have dbi/dXj = —1, so that 

The coordinates (Xj, Zj) of the moving points turn out to satisfy the Hamiltonian system 

■ dH dH 

X, - z 3 - — . (5.15) 
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Strictly speaking, it is Xj and log Zj that play the role of Darboux coordinates. 

The algebraic relations (j5.12j) imply that this Hamiltonian system is separated to a 
collection of independent systems on a two-dimensional phase space with coordinates 
(Xj,Zj); this is the central concept in "separation of variables" [IBIIHEE]- An integrable 
Hamiltonian system with a maximal number of constants of motion Ii, ■ • ■ , I n in involution 
is separable if it has a suitable set of Darboux coordinates Xj,fij that satisfy a set of 
equations of the form 

fj(\j, [Mj, h, ■ ■ ■ , I n ) = (j = 1, - • • , n). 

Geometrically, these equations define a family of Liouville tori, which are thus "separated" 
to a direct product of one-dimensional level curves. For most integrable Hamiltonian 
systems with a Lax pair, these equations take an identical form 

f(\j,fij, h,---,I n ) =0 (j = 1, -,n), 
which stems from the the equation 

f(X,fi,h,---,I n ) =0 
of an associated spectral curve, as we have just illustrated for periodic Toda chains. 



5.2.2 Periodic dressing chains 

Bearing in mind this description of periodic Toda chains, we now turn to periodic dressing 
chains. The algebraic equations (|5.7jl can be solved for I/s, I = 1, • • • , g, as 

» Zj + Q(X j )z~ 1 - 7 Af db t 



B'(Xj) dX j 



for N = 2g + 1, and 



^ z J + Q(X J )z- 1 -2Xf 1 -I X° db i 
U = "S »M 9X~ (5 ' 17) 

for = 2g+2. Each Ii thus becomes a function of the 2g variables Ai, • • • , X g and Zi,---,z g 
(as well as of x, through the x-dependence of bo and Io, provided that a ^0). In view of 
the case of periodic Toda chains, it seems likely that the first nontrivial coefficient 1\ will 
play the role of the Hamiltonian. We have indeed the following result, which says that the 
analogy is essentially valid up to a simple factor. Note, however, that the Hamiltonian 
system thus obtained is non- autonomous (and presumably not separable) in the case where 
a ^ 0. 
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Proposition 2 \5.§) and A5.9\) can be cast into the Hamiltonian syst 



em 



where 



bo 

More explicitly, the Hamiltonian takes the form 



H=^-. (5.19) 



for N = 2g + 1 and 



r ^ , J+ g(A J )^-2Af 1 -/ Af (52i) 



forN = 2g + 2. 

Proof. It is obvious from the explicit form of H that 

dH _ Zj - QiX^zJ 1 



dzj B'(Xj 



The first equation of ()5.18|) thus turns out to be the same thing as ()5.8|) . To derive the 
second equation of (|5.18|) . we notice that the equations 

z 2 k - P(X k )z k + Q(X k ) = 

are identically satisfied if U are understood to be a function of A/s and z/s defined by 
(|5.16j) or (|5.17j) . Differentiating these equations by Xj yields the equations 

{-P\X k )z k + Q'(X k ))8 jk - -*fK~ l *k = o. 

t=i aA o 

We can solve these equations for dlj/dXj by the inversion formula (|A.5|) as 

dh ^ (-P'(X j )z k + Q'(X j ))6 jk db t 



E 



dXj j^J B'(X k )z k dX k 

P\X 3 )z 3 -Q\X 3 ) db £ 
B'(X 3 )z 3 dX/ 
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Since bi = — 6 (Ai + • ■ • + X g ), we have 

dh 

so that 

dH Zj dh _ P'iX^zj - Q'(Xj 



' J d\j~ b dXj B'(\j) 



The second equation of (|5.18|) thus turns out to take the form 

P'{X 3 )z 3 - Q'(Xj) 



Z 3 



B'(Xj) 

which is exactly ()5.9|) . □ 



5.3 Reconstructing Lax equation 

A significant feature of the spectral Hamiltonian system (j5.18j) is that it contains no infor- 
mation on P(X) (except for I ). In other words, this Hamiltonian system has apparently 
"forgotten" what the spectral curve is. The only external input of the Hamiltonian are 
the polynomial Q(X) and the two quantities 60 and Iq that are a constant or a simple 
combination of v. 

Actually, alongside the spectral curve, one can reconstruct the transition matrix T(A) 
itself from the Hamiltonian system. Thus the spectral Hamiltonian system is actually 
equivalent to the Lax equation ()4.4|) . 



5.3.1 If iV is odd 

We first consider the case where iV = 2g + 1. Since bo = 1, the external input of the 
Hamiltonian are Q(X) and Io = 2t>; v is assumed to be a function that satisfies ()2.6j) . 

The reconstruction proceeds as follows. 

i) Since b = 1, we set 

B(X) = X 9 + hX 9 - 1 + ■ ■ ■ + b g = IJ(A - Xj), (5.22) 

and define the polynomial 

A(X) = a X 9 + aiX 9 ' 1 + ■ • • + a g , a = ^, (5.23) 
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by the equations 



M*i) = *j (; = !,...,»). (5-24) 

ii) We now define 

D(X) = P{X)-A(X), (5.25) 

C(A) = mm-m (5 26) 

D(X) is a polynomial of the form D(X) = d^X 9 + dxX 9 " 1 + ■ ■ - + d g with do = Io — a = Iq/2. 
C(X), too, becomes a polynomial, because the numerator vanishes at the zeroes A = Xj 
of the denominator, see (|5.3|) . Since A(A)-D(A) — Q(A) is a polynomial of degree 2g + 1 
with the leading coefficients being equal to 1, C(A) is a polynomial of the form C(X) = 
X 9+1 + C\X 9 + • • ■ + c 9+ i. The four matrix elements of T(A) have been thus determined. 

iii) Let us examine the first equation of ([5. 18)1 . Calculating the right hand side explic- 
itly, we can rewrite this equation as 

B'(X)Xj = A(Xj) - D(Xj). 

By the identity B(Xj) + B'(Xj)Xj = 0, the last equation can be further rewritten as 

B(Xj) + A(Xj) - D(Xj) = 0. 

On the other hand, since bo — 1 and a = do, B(X) + A(X) — D(X) is a polynomial of degree 
less than g. Such a polynomial cannot vanish at g distinct points unless it is identically 
zero. We thus arrive at the equation 

B(X) + A(X) - D(X) = 0, 

which is exactly (J4.6|) . 

iv) Let us examine the second equation of ()5.18j) . Viewing the algebraic relations 
z\ — P(Xk)zk + Q(Xk) = as identities satisfied by J/s, we repeat the calculations in the 
proof of Proposition El and obtain the expression 

OH _ P'iX^Zj - Q'(Xj) 

Zn 



3 dXj B'(Xj) 

for the right hand side of the second equation of (|5.18|) . On the other hand, differentiating 
the same identities by x along the trajectories of (|5.18|) gives 

(2 Zj - P'{X 3 )) Zj + {-P'iX^Zj + g'(A,))A J - PiX^Zj = 0, 



22 



which implies 



P'(A^ - Q'(X,) • P(A,) 
2z i -P(A i ) Ai+ 2,,-P(A i ) % 



P'(A,)^ - Q'(Xj) P(A 



//'(A.;) 2z j -P(X j ) 3 

For the second equation of (|5.18j) to hold, the algebraic relations P(Xj) — (j = 1, • • • , g) 
have to be satisfied This implies that P(A) is divisible by B(X). Moreover, since P(A) is a 
polynomial of degree g with the leading coefficient being equal to a, we eventually obtain 
the equation 

B(X) a ' 

which is nothing but ()4.19|) . 

v) In order to derive f)4.5|) . we start from the identity A(Xj) + A'(Xj)Xj = Zj that can 
be obtained from ()5.24j) . Since Xj and z 3 - can be calculated as shown in (|5.2|) and 
we have 

■ D'(Xj)A(Xj) + A'(Xj)D(Xj) - Q'(Xj) 

{ j) ~ B\X~) ' 

On the other hand, the value of C(A) at A = Xj can be calculated by de L'Hopital's 
formula as 

c ( A,) = lim *m?tzm 

v 3> B(X) 

A'(X j )D(X j ) + A(X j )D'(X j )-Q'(X j ) 
B\X 3 ) 

Since the right hand side of these two equalities are identical, we find that 

A{X 3 ) = C(Xj) 

for j = 1, • • •, g, so that A(X) — C(A) is divisible by -B(A). On account of the Laurent 
expansion of the quotient at A = oo, we can deduce that 

A(X)-C(X) x . 

r> / w = -A + a - Cl + 6i. 

This equation coincides with ()4.5|) upon identifying 

ui = —d + Ci — b\. (5-27) 
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vi) The final step is to verify ()4.7|) and (J4.8|) . Firstly, we now have ()4.5|) and (|4.19j) at 
hand, and using these equations we can derive (|4.8|) : 

D(X) = P(A) - A(X) = (A + Ul + a)B{\) - C{\). 

To derive (|4~7|l . we differentiate the identity Q(A) = A(A)D(A) - B(X)C(X) by x and 
eliminate the derivatives A(X), B(X) and D(X) using the differential equations that have 
been proven: 

= A(X)D(X) + A(X)D(X)-B(X)C(X)-B(X)C(X) 
= B(X)((X + a + Ul )A(X) - (A + m)D(A) - (7(A)). 

Removing the common factor B(X), we obtain (|4.7jl . 



5.3.2 If iV is even 

Let us now consider the case where = 2g + 2. The external input of the Hamiltonian 
are Q(X), b = v and I = v 2 + • • • + fiig+i] v is assumed to be a function that satisfies 

flU). 

The reconstruction of T(A) is mostly parallel to the previous case. -B(A) is the poly- 
nomial 

B(X) = b X 9 + hX 9 ' 1 + --- + b g = b f[(X- A,). (5.28) 

A(X) is a polynomial of the form 

A(X) = X 9+1 + a X 9 + ■ ■ ■ + a g , a = \{h ~ &o), (5.29) 
and defined by the equations 

A(X j ) = z j (j = l,---,g). (5.30) 

C(X) and -D(A) are given by 

D(X) = P(X)-A(X), (5.31) 

C W = (5,2) 

They become polynomials of the form -D(A) = A 9+1 + c?oA 9 + d\X 9 ~ l + ■ ■ ■ + d g and 
C(A) = cqA s+1 + C\X 9 + • — h c 9+ i with <i = (-^o + &o)/2 and Co = &o- As we have done for 



24 



odd iV's, these four polynomials turn out to satisfy (|4.5|) . (|4.6|) . (|4.7|) and (|4.8|) with iti of 
the form 

-a + ci + 6i 

Wi = . (5.33) 



We have thus proven the following: 

Proposition 3 Given a solution of the spectral Hamiltonian system A5.1ty) . one can re- 
construct the transition matrix T(A) and the potential u\ by A5.2fy) - \5.2l\) (for N = 
2g + l) and UTM) - JO^) (for N = 2g + 2). 

The contents of this section can be summarized as follows. 

Theorem 1 The Lax equation \4-4\) °f the transition matrix T(A) is equivalent to the 
spectral Hamiltonian system A5.18\) in the spectral Darboux coordinates. 



6 Inverse problem and Poisson structures 

We now turn to the "inverse problem" , namely, the problem of reconstructing the phase 
space coordinates of the periodic dressing chain, or of the Noumi-Yamada system, from 
the spectral Darboux coordinates. We can solve this problem and find that the spectral 
Hamiltonian system and the Noumi-Yamada system are connected by a locally invertible 
rational map (Ai, . . . , X g , zi, . . . , z g , v) i— > (fx, ... , f N ). Moreover, this map turns out to 
be a Poisson map connecting two odd- dimensional Poisson structures. 



6.1 General setup and case study 

The phase space of the spectral Hamiltonian system has a natural symplectic (hence, 
even- dimensional Poisson) structure. Actually, one can extend this phase space to an 
odd-dimensional Poisson manifold with an extra dimension. On the other hand, the 
phase space of the Noumi-Yamada system is known to have an odd-dimensional Poisson 
structure [UHE]. 

The two Poisson structures are formulated as follows. 

1. The first Poisson structure is defined on the extended phase space of the spectral 
Hamiltonian system (j5.18)l with the auxiliary variable v added as an extra dimension. 
This is a (2g + l)-dimensional Poisson manifold with the Poisson brackets 

{Xj, z k ] = 5 jk , {Xj, X k } = {zj,z k } = {Xj,v} = {zj,v} = 0. (6.1) 
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In particular, v is a Casimir function, i.e., a central element of the Poisson algebra 
of coordinates. 

2. The second Poisson structure is defined on the phase space of the Noumi-Yamada 
system endowed with the Poisson brackets 

{fmy fn} $m+\,n ^m—l,n- (6-2) 

Since the variables /„ are required to satisfy the constraint (|2.18|) , the phase space 
in this case, too, is (2g + l)-dimensional. The auxiliary function v — (/i + - ■ -+/jv)/2 
is a Casimir function of this Poisson structure as well. 

Let us now examine the simplest two cases, i.e, N = 3 and N — 4. 
6.1.1 Case study: N = 3 

The three-periodic dressing chain has the three fundamental variables V\, v 2 , v 3 , which are 
connected with the variables of the Noumi-Yamada system of A 2 type: 

/i = vi + v 2 , f 2 = v 2 + v 3 , f 3 = v 3 + vi. (6.3) 

The phase space of the Noumi-Yamada system is now a three-dimensional Poisson mani- 
fold with the Poisson brackets 

{/l,/ 2 } = {/2,/8} = {/3,/l} = -l (6.4) 

or, equivalently, 

{v h v 2 } = -1, {vi,v 3 } = 1, {^2,^3} = -I- (6.5) 

It is easy to verify that the auxiliary function v = v\ + v 2 + v 3 = (fi + f 2 + f 3 )/ is indeed 
a Casimir function of this Poisson structure. This is the Poisson structure that appears, 
along with a generalization to odd N's, in the work of Veselov and Shabat [26J. 
The equation of motion (j2.14|) can be rewritten as 

A = {fu H}, h = {f 2 , H}, f 3 = {f 3 , H} + a, (6.6) 

where H is the function that can be calculated from the transition matrix as 

H = Tr T(0) = hhh + Pih + Mi- (6-7) 
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Actually, this is the same Hamiltonian as that of (|5.18|) . i.e., H = I\. This "Liouville 
form" of the equation of motion is derived from a quite different point of view in the work 
of Noumi and Yamada ^Hj 

The construction of the spectral Darboux coordinates now takes a particularly simple 
form, because B{X) and A(X) are a linear function of A: 

B(X) = A + f x f 2 + ft, A(X) =vX + (A/ 2 + ftK (6.8) 

The spectral Darboux coordinates Xi,zi are thus determined explicitly as 

Ai = -/i/ 2 -ft, Z! = hXi. (6.9) 
Note that this expression of Ai and z\ can be solved for f\ and f 2 as 

, z\ , Ai(Ai + 

J2 = t-, Ji = , (6.10) 

so that, along with f 3 = v — fi — /2, all / n 's are a rational function of X\,z\ and v. 

Moreover, one can easily verify that the Poisson brackets of are consistent with the 
Poisson brackets 

{X 1 ,z 1 } = z 1 , {X l ,v} = {z 1 ,v} = (6.11) 

of Ai, Z\ and v. In other words, the rational map (fi, f 2 , ^3) 1— > (Ai, Zi,v) is Poisson. 
As a final remark, let us point out that A(X) can be also expressed as 

A(X) = B(X)vi + f 2 X. (6.12) 

We shall encounter a similar expression in the case study for N = A below. 



6.1.2 Case study: N = 4 

The four-periodic dressing chain has the four fundamental variables V\, . . . , V4 under the 
quadratic constraint (jB.l|) . The four variables fx , . . . , f± of the Noumi- Yamada system of 
A<P type are defined by the linear relations 

fx = Vi + v 2 , f2=v 2 + v 3 , f 3 = v 3 + t> 4 , fi = v A + vi, (6.13) 

and obey the constraint fi + f 3 = fi + fi- They are thus a set of (redundant) coordinates 
of a three-dimensional Poisson manifold with the Poisson brackets 

{h, h} = {/ a , h] = {/ 3 , h} = {U, h} = -1, 

{/i,/s} = {/i,/4} = {/ a ,/4} = 0. (6.14) 
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The auxiliary variable v = fx + f 3 = f 2 + f& is a Casimir function. 

As we have seen in Section 2, the linear relations connecting / n 's with t> n 's can be 
solved for u n 's under the constraint (jB.l|l . The outcome is a rational expression of t> n 's in 
terms of /„'s, such as 

vi = ^-(2v 2 - 4/2/3 + oci - a 2 + a 3 - a 4 ), (6.15) 

etc. The Poisson brackets of t> n 's are thereby determined by the Poisson brackets of / n 's 
as 

r 1 fs r -1 fl . /} r -1 fl 

{Vx,V 2 \ = , {vx,v 3 } = + — , {vx,v 4 } = — , 

V V V V 

{^2,^3} = - — , {v2,v 4 } = -— + — , {^3,^4} = - — • (6.16) 

V V V V 

For example, the Poisson bracket of v 1 and v 2 can be correctly derived from the identity 
{vx, v 2 } = {vx, fi - vx} = {vx, fx} and (|H.lftjl. 

The equation of motion ()2.17|) again turns out to take Liouville form 

ft = {fi, H}, f 2 = {/ 2 , H}, h = {/a, H}, U = {U, H} + a (6.17) 
with the Hamiltonian 

H = v- 1 Tr T(0) = v-\f x f 2 hU + A/3/4 + A/1/4 + A/1/2 + PM- (6.18) 

This Hamiltonian H is the same thing that appears in (|5.18|) . i.e., H = Ix/v. 

The construction of the spectral Darboux coordinates is mostly parallel to the case of 
N = 3, because B(X) is a linear function of A. More precisely, we have 

5(A) = UA + /J2/3+ A/3 + &/i, 

A(X) = B{\)vx + (/ 2 / 3 + A + f3 2 )\. (6.19) 

Note that we have expressed A (A) in a particular form, which we shall generalize to all 
iV's. The spectral Darboux coordinates Ai,£i are thus determined as 

Ai = -v-\hf 2 h + A/s + A/i), = C/2/3 + Ax + A)Ai- (6.20) 

Combined with the linear relations fx + f 3 = f 2 + f 4 = v, these formulae of \x,zx can 
be solved for / n 's. This yields an expression of / n 's as a rational function of Xx,zx and 
v, which can be assembled to a rational map (Ai, zx,v) 1— > (/1, / 2 , f 3 , fi). This map, too, 
turns out to be a Poisson map with regard to the aforementioned Poisson structures. 
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6.2 Solving inverse problem 

Let us now consider, in general, how to reconstruct t> n 's or / n 's from the spectral Darboux 
coordinates Xj, Zj and the auxiliary variable v. Since the latter carry the same information 
as the transfer matrix T(A), the problem is to reconstruct t> n 's or f n 's from the transfer 
matrix. 

To solve this problem, we make use of a new auxiliary polynomial A(X) such that 

A(X) = B(X)vi + A(X)X. (6.21) 

This polynomial and B(X) have several interesting properties, which eventually lead to a 
solution of the inverse problem. 

6.2.1 Structure of A(X) and B(X) 

The existence of the polynomial A(X) has been hinted in the case study for N = 3 and 
N = 4. Not only the existence for all N > 3, we can also deduce the following result on 
its precise structure. Note, in particular, that A(X) and B(X) turn out to be a function 
of / n 's rather than t> n 's. 

Proposition 4 For all N > 3, there is a polynomial A(X) that satisfies 116. 21)) . Moreover, 
this polynomial can be written as 

N - 2 ( d 2 \ 

MX) = II ! + + rlf (/2-/ ff -i) (6-22) 

n=2 V OjnOjn+lJ 

for N > 4 and A(X) = f% for N = 3. Moreover, B(X), too, can be similarly expressed as 

N-2 / Q2 \ 

B{\) = II 1 + ( A + ^) ^7 (A-U ( 6 -23) 

n=l V OJnOJn+lJ 

Proof. Let us note that A (A) and -B(A) can be expressed as 
A(A) 



Tr( 1 0)V N {X + p N - 1 )---Vi{X) ( J j 
Tr(l 0)v A r(A + /3 JV _ 1 )---F 2 (A + /3 1 ) f 



X + v\ 



and 



5(A) = Tr(l o)v N (X + f3 N _ l )---V 1 (X)(^° i 

= Tr(l ) V N (X + P N ^ 1 )---V 2 (X + (3 1 ) f 1 j 
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Substituting 



we can rewrite the last expression of A(X) as 



A(X) = Tr(l o)V N (X + (3 N „ 1 )---V 2 {X + (3 1 )(^ 



Vl 



Tr( 1 ) V N (X + Pn-i) ■ ■ ■ V 3 (X + fo) | J ) A. 

Since the first term on the right hand side is equal to B(X)vi, one obtains 1)6.21)1 upon 
defining A (A) as 

1 

V-2 



A(A)=Tr 1 o \ Vn (x + (3 n _ 1 )...V 3 (X + P 2 ) 



We can now apply the method used for evaluating the trace formula ()4.23|) to the right 
hand side of these formulae. This yields (j6.22j) and (|6.23J) . □ 

6.2.2 Three-term recursion relations 

The foregoing expression of A(X) and B(X) shows that they are the first two members of 
the finite sequence of polynomials 

N-2 / Q2 \ 

F m (X)= II ! + +fl0-~77w— )(fm---f N -i) (m = l,...,iV-2) (6.24) 

n=m \ OfnOfn+l/ 

and i 7 jy_ 1 (A) = Jn-i- An interesting property of these polynomials is that they satisfy a 
set of three-term recursion relations as we show below. It deserves to be mentioned that 
a similar finite sequence of polynomials take place in Moser's method for solving finite 
non-periodic Toda chains [TT] . 

Proposition 5 The polynomials F m (X) satisfy the three-term linear recursion relations 

F m (X) = f m F m+1 (X) + (A + (3 m )F m+2 (X) (6.25) 
and the auxiliary relation 

= F m+1 (X). (6.26) 
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Proof. The definition of F m (X) can be rewritten as 



F m (\) 



1 + (A + /3 m ) 



/m-^m+1 

(X) + (X + /3 m ) 



dfmdfm+l J 

dF m+1 (X) 



(fmF m +l(X)) 



m+l 



It is also obvious from the construction that F m+ i(X) does not contain f m . Since f m 
appears only in the first term on the right hand side, one immediately finds that ()6.26j) 
holds. The last equality thereby turns into the three-term linear recursion relation. □ 

An immediate consequence of the three-term recursion relations is that F m (X) has the 
determinant formula 



F m (X) 



fm 
-X-Pr, 



1 

/m+l 



-A-/3 



7V-2 



1 

/jv-i 



(6.27) 



This formula shows that the polynomials F m (X) are also closely related to the Lax pair of 
Noumi and Yamada. It is easy to prove this formula: Firstly, one can easily verify that 
these determinants satisfy the same recursion relations. Secondly, the formula holds for 
m = N — 2 and m = N — 1. Consequently, the determinant has to coincide with F m (X) 
for all m's. 



6.2.3 Reconstructing / n 's 

We are now in a position to solve the problem of reconstructing / n 's from the spectral 
Darboux coordinates Xj, Zj and the auxiliary variable v. A clue is the relation 

F m (—(3 m ) = fm,F m +i(—Pm) (6.28) 

that can be derived from the three-term recursion relations. This relation may be thought 
of as a formula that determines f m from F m (X) and F m+1 (A). Similarly, -F m+2 (A) can be 
expressed as 

7-i /\\ F m (X) — fmF m +l(\) , . 

i t m+2\A) — : n ■ [b.AU) 

A + Pm 

Thus f m and -F m+2 (A) can be recursively determined from Fi(X) = B(X) and F 2 (X) = 
A(X). As regards B(X) and A(X), they are polynomials of the form 

B(X) = b X 9 + --- + b g , A(X) = a X 9 + •■■ + %, 
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where b g = 1, ao = for N = 2g + 1 and bo = v, ao = 1 for N = 2g + 2, and satisfy the 
algebraic relations 

B(Xj) = 0, ^ = A(Xj)Xj (6.30) 
as one can readily see from ()6.21|) . Consequently, by the interpolation formula ()A.2|) . 



B(A)=fcn(A-A,), ^)=(| + g B , (Aj)( ; j _ Aj)Aj . )^(A). (6.31) 
Starting from these polynomials, one can thus determine / n 's successively as 

/l = 2tM, /2= + etc. (6.32) 

^(-A) B(-A) - JU(-ft) 

In particular, each / n is a rational function of the spectral Darboux coordinates and t> . 

This completes the reconstruction of / n 's in the case where N = 2g + 1. If iV = 
2g + 2, we have to confirm that / n 's satisfy the constraint f)2.18j) as well. Actually, this 
is an immediate consequence of the recursive formulae of F m (A)'s: Starting from the 
aforementioned form of -Fi(A) = B(X) and i^(A) = ^4(A), one can show by induction that 

-^2n(A) = X g ~ n+1 + terms of lower degree, 
^2n+i(A) = (v-fi-fa /2n-i)A 9 ~ n + terms of lower degree (6.33) 

for n = 1, 2, . . . , g, so that, eventually, 

/as+i = F 2g+ i(X) = v - /i - f 3 - . . . - f 2g -i- 



This implies fl2T8j) . 

We have thus proven the following. 

Theorem 2 The Noumi- Yamada system is connected with the spectral Harailtonian sys- 
tem by a locally invertible rational map (Ai, . . . , X g , z±, . . . , z g , v) h-> (f\, . . . , /at). This 
map consists of rational functions recursively calculated as 116. 



6.3 Comparing two Poisson structure 

Now that the spectral Hamiltonian system and the Noumi- Yamada system are con- 
nected by a locally invertible map, one can compare the Poisson structures on the (ex- 
tended) phase spaces of these systems. Actually, it is more convenient to consider this 
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issue at a place in between, namely, on the (2g + l)-dimensional space with coordinates 



If one reviews the foregoing construction carefully, one will soon notice that the map 
constructed above can be factorized to the composition of two maps as 

(Ai, . . . , X g , z x , . . . , z g , v) A (ai, . . . , a g , h, . . . , b g , v) A (f t , . . . , f N ), (6.34) 

each map being defined by (|6.31|) and (|6.32|) respectively. Moreover, whereas 71 a locally 
invertible polynomial map, 72 is birational; 72 -1 is a polynomial map defined by (|6.22j) . 



The Poisson structures (jfi.lj) and ()6.2|) induce, via these maps, two Poisson structures 
on the space in the middle. We show bellow that the two induced Poisson structures are 
the same. This implies that the Poisson structures on the both ends are locally equivalent. 

6.3.1 Poisson algebra induced by (16. lj) 

We first examine the Poisson structure induced by (|bM|) . Since v is obviously a Casimir 
function, the problem is to calculate the Poisson brackets of a/s and b/s. This calculation 
can be done most neatly in terms of generating functions, i.e., A(X) and B(X). The 
outcome is a kind of "quadratic Poisson algebra" as follows. 

Proposition 6 The Poisson structure of \6.1\) induces a quadratic Poisson algebra of 
A(X) and B{X). The induced Poisson brackets read 



&i, . . . , a g , bi, 




§M and v = (A + • • • + f N )/2. 



{A(X),A^)} = {B(X),B^)} = 



A — [L 



(6.35) 



where A and fi are understood to be independent parameters. 



Proof. The proof consists of several steps. 

i) Since {A^, Xk} = {Xj, v} = 0, it is obvious that {B(X), -B(//)} = 0. 

ii) Let us substitute Zk = A(Xj)Xj in the fundamental relations {Xj, Zk} = SjkZk- Since 



{Xj, A(X k )} = {Xj, A(X)}\\=x k + {Xj, X k }A'(Xk) = {X v A(X)}\ X=K 



this yields the equalities 



{Xj,A(X)}\ x=X: 



■k 



SjkA(Xk) 
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for k = 1, . . . , g. The interpolation formula (jA.l|) can be used to assemble them to a single 
relation: 



{Xj, A(X)} 



i(A,)5(A) 



B'iXj^X-XjY 

iii) We now consider the obvious equality {B(Xj), A(X)} = 0. The left hand side can 
be expanded as 

{B(X j ),A(X)} = {£(//), I(A)}| M=A . + {Xj,A(X)}B'(Xj). 

Therefore, applying the last formula of step ii to the second term on the right hand side, 
we obtain the equalities 

M^)B(X) 



{B(fi),A(X)}\, =Xj = -{A„A(A,)} J B'(A J ) 



A-A, 



for j = 1, . . . , g. Now the interpolation formula is invoked once again to pack them into 
the single relation: 

{b M ,a(x) } - -E B , (A . )(A _ A . )(M _ A . r 

Substituting 



(A — Xj)(jj, — Xj) X — /i — Xj X — Xj J ' 
we can split the right hand side of the last equality into the difference of 

* A(Xj) _ A(X) d 
j^x B'(Xj)(X - Xj) B(X) b 

and the same quantity with A replaced with \i. This yields the expression of {v4(A), B(fi)} 
in JOSJ). 

iv) Let us now rewrite the fundamental relations {zj,z k } = in the same way. Sub- 
stituting Zj = A(Xj)Xj and z k = A(X k )X k , we encounter the quantity 

{A{Xj),A(X k )} = {X j ,A(X k )}A'{X j ) + {A{X),X k }\ x=Xj A'{X k ) 
+ {A(X),A(»)}\ 

\=\j ,[J,=\k • 

The first term on the right hand side is equal to Sj k A(X k )A'(X k ), which cancels with the 
second term, because the latter can be calculated as 

{A(X),X k }\ x=Xj A'(X k ) = - hm ^^£^Li'(A fc ) = -6 jk A(X k )A\X k ). 

We can thus deduce that {A(X), A(fi)}\ x=Xj ^ =Xk — for j, k = 1, . . . , g. This implies that 
{A(X),A(ti)} = 0. □ 
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6.3.2 Poisson algebra induced by (J6.2|) 

We now calculate the Poisson brackets of A(A) and B(X) with regard to the Poisson 
structure of (|fi.2j) . The three-term recursion relations play a key role here again. 

Proposition 7 The polynomials F n (X), n = 1, ... ,N — 2, pairwise satisfy the Poisson 
relations 

{F n+1 (A),F„ +1 (/i)} = {F n+2 (X),F n+2 (fi)} = 0, 

si? (\\ t? (,\\ F n+2WF n+1 {n) - F n+1 (X)F n+2 (fi) . , 

{F n+ i(A),F n+2 (/i)} = (6.36) 

wi^/i regard to the Poisson structure of \b\ §) . In particular, this induces the same Poisson 
algebra \6. Sty as the one induced by the Poisson structure of \6.1)) . 

Proof. We prove ()6.36|) by induction that starts from n = N — 3 and proceeds downward 
to n = 0. In the case where n = N — 3, ()6.36|) is a statement for -F/v_ 2 (A) = /jv_2/jv-i + 
A + Pn-2 and Fjy-i(A) = /jv-i) which one can directly confirm. Let us now assume ()6.36|) 
and show that the same Poisson relations with n replaced by n — 1 are satisfied. More 
explicitly, our task is to derive the Poisson relations 

rp / n T? ( W n (T? f\\ T? I W Fn+lWF n (fi) - F n (X) F n+1 (fl) 

{F n {\), F n (fi)\ = 0, {F n (X),F n+1 (fi)\ = (6.37) 

A — n 

of the next stage from (|6.36j) . This is indeed achieved by substituting 

F n (X) = f n F n+1 (X) + (A + (3 n )F n+2 {X), F n (n) = f n F n +M + + (3 n )F n+2 {fi), 

and doing calculations. Let us illustrate this procedure for the second part of ()6.37j) . 
Upon substitution, the left hand side of this Poisson relation can be expanded as 

{F n (A),F n+1 (/i)} = {f n ,F n+1 (fx)}F n+1 (X) + f n {F n+1 (X),F n+1 (fi)} 

+ (A + /3 n ){F n+2 (A),F n+1 (/i)}. 

The first Poisson bracket on the right hand side can be calculated as 

{fn,F n+ M} = - d -^±M_ = - Fn+ M. 

The other Poisson brackets can be handled by the Poisson relations (J6.36|) that we have 
assumed. One can now confirm the first part of (J6.37|) by straightforward calculations. 
The first part of ()6.37|) can be verified in much the same way. □ 

We have thus arrived at the following result. 
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Theorem 3 The rational map connecting the spectral Hamiltonian system and the Noumi- 
Yamada system is a Poisson map with regard to the Poisson structures defined by \6.1\) 
and \6. ty) . 

7 Conclusion 

Although Adler's Lax pair for periodic dressing chains has an unusual structure, we have 
been able to apply the usual method for isospectral and isomonodromic deformations with 
minimal modifications. Firstly, the spectral curve and the spectral Darboux coordinates 
can be constructed from the transition matrix around the periodic chain. Secondly, the 
Lax equation of the dressing chain can be converted to a Hamiltonian system — the 
spectral Hamiltonian system — in the spectral Darboux coordinates. Thirdly, the dy- 
namical variables of the dressing chain can be reconstructed from the spectral Darboux 
coordinates by an algebraic procedure. As a byproduct, we have been able to confirm 
that the natural symplectic structure of the spectral Hamiltonian system is consistent 
with one of the previously known Poisson structures of the periodic dressing chain or of 
the Noumi-Yamada system. 

It is significant that our "spectral" description of the Hamiltonian structure has turned 
out to be closely connected with Noumi and Yamada's Lax pair. This seems to indicate 
that one can reformulate everything in the language of Noumi and Yamada's Lax pair. 
In this respect, a very intriguing problem is to develop a similar spectral description of 
the sixth Painleve equation using Noumi and Yamada's new Lax pair |2Uj . 

It should be stressed that we have been able to digest just a small part of the full 
list of various dressing or Darboux chains. For instance, Adler's work [2] deals with some 
other Painleve equations. Moreover, Recent work of Willox and Hietarinta |2Z| seems to 
cover a broader class of Darboux chains. 
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A Interpolation formulae 



Let B(X) be a polynomial of the form 

B(X)=j2btX g - t = b f[(X-X s ), 

and assume that the zeroes Xj are mutually distinct. The simplest form of Lagrange's 
interpolation formula reads 

F(X) ° F(Xj) 

B(X) ~t B'(Xj)(X - Xj) ' 1 • } 

which holds for for any polynomial F{X) of degree less than g. If the degree of F(X) is 
greater than or equal to g, extra terms appear on the right hand side as 

where G(X) is a polynomial. 

Let us consider, in particular, the case where F(X) — X e , £ — 0, • • • , g — 1. This leads 
to the identities 

B(X) ~t B'(Xj)(X - Xj) 
Picking out the residue of both hand sides at A = oo gives the formula 

which has been used in the theory of finite-band integration. 

Another application is an inversion formula of a system of linear relations with Van- 
dermonde coefficients. Consider the linear relations 

E&Af = ifc U = l,-~, 9)- (A.4) 

The problem is to solve these equations for £/s. A clue is the auxiliary function {(A) = 
Z)?=i6A 9_£ - This is a polynomial of degree less than g, and satisfies the interpolation 
relations £(Xj) = rjj for j — 1, • • • , g. The the interpolation formula thereby yields 

ro.) V ^ B{X) V ^ dBW 

^ > ~1 B'(Xj)(X - Xj) - UB'(Xj) dXj ■ 

Extracting the coefficients of the polynomials on both hand sides gives the inversion 
formula 

6 = -g5feS ( ' = 1, - ,9) - <A - 5) 
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B Equivalence of even-periodic dressing chain and 
Noumi-Yamada system 

The identity (|2.15jh combined with (|2.3J1 . yields the quadratic constraint 

2g+2 

E (- 1 ) n (^n+i -vl + an) = Q (B.l) 

n=l 

on VnS. This constraint can be rewritten as 

g+l 1 2g+2 

fc=l Z n=l 

Substituting f| fc — fffc-i = ( u 2fc — u 2fc-i)/2k-i> the quadratic constraint turns into the linear 
constraint 

g+l 2g+2 

EK - U2fc-l)/2*-l = ^ E ("l) n «n- (B.2) 
fc=l ^ n=l 

One can use this linear constraint to remove the ambiguity for solving (j2.12j) as follows. 
Firstly, solve (j2.12|) for v 2 , v 3 , . . ., successively as 

V2 = fi~ v x , v 3 = f 2 - v 2 = f 2 - fi + «i, . . . , 

V n = fn + l-fn + --- + (-1)" fl " (-1)^1- (B.3) 



One can thereby eliminate v^s from (|B.2)1 and obtain the equality 

(g+l \ 2 g+l 1 2 S +2 

E/«-i " 2 E /2i/ 2fc +i-2^t;i/ 2fc _ 1 = -E(-ira n . 

\fe=l / l<i<fc<9 fc=l «=1 

Solving for v\ gives the expression 

1 / 1 29+2 \ 

v, = -U-2 J2 Mu+i ~ o E ("!)"«» (B.4) 

ZU V l<J<fe<9 Z n=l / 

of t>i as a rational function of / n 's and v. This, in turn, determines the other v n 's by 
(JESJ). 

Thus, although the linear relation (J2.12j) itself cannot be solved for t> n 's uniquely, 
the quadratic constraint (|B.1J) removes the ambiguity and yields a rational expression of 
Vi, . . . , v 2g +2 in terms of / n 's and v. In other words, the (2g + 2)-periodic dressing chain 
is connected with the Noumi-Yamada system of ^29+1 tyP e by a birational change of 
coordinates. 
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